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ABSTRACT 


This  report  summarizes  the  work  done  by  the  SDL  during  the 
period  April  through  June  1968,  and  is  primarily  concerned  with 
seismic  research  activities  related  to  the  detection  and  identi 
fication  of  nuclear  explosions  and  earthquake  phenomenon.  Also 
discussed  are  the  support  tasks  and  data  services  performed  for 
other  participants  in  the  VELA-Uniform  project. 


I .  INTRODUCTION 


This  quarterly  report  summarizes  the  technical  work,  sup¬ 
port  effort,  and  service  tasks  completed  during  the  period 
April  through  June  1^68.  Current  or  past  work  is  mentioned 
only  if  it  relates  to  the  present  discussions. 

Reviews  of  technical  reports  completed  during  the  report¬ 
ing  period  are  contained  in  Section  II  under  descriptive  head¬ 
ings.  Section  III  is  a  summary  of  the  support  and  service 
tasks  performed  for  in-house  projects  and  for  other  VELA- 
Uniform  participants.  This  report  concludes  with  Appendix  A 
in  which  are  listed  the  organizations  who  received  SDL  data 
services  during  the  period. 

II.  WORK  COMPLETED 

A.  Precision  Location  of  Underground  Nuclear  Explosions 
Using  Teleseismic  Networks  and  Predetermined  Travel- 
Time  Anomalies 

This  report  is  primarily  concerned  with  the  teleseismic 
location  accuracy  obtainable  for  a  particular  region  by  using 
either  station  residuals,  travel-time  anomalies,  or  no  correc¬ 
tions  at  all.  Secondarily,  we  are  concerned  with  the  perform¬ 
ance  of  two  computer  programs,  LOCATE  and  SHIFT,  presently  in 
use  at  the  Seismic  Data  Laboratory.  The  principal  difference 
between  the  two  programs  is  that  SHIFT  minimizes  in  a  least- 
squares  sense  the  relative-anomaly  errors  rather  than  the  ab¬ 
solute  residuals. 

Also  investigated  is  a  technique  whereby  limits  of 
the  solution  are  estimated  either  on  the  basis  of  an  acceptable 
network  standard  deviation  of  time  errors  or  of  maximum  rela¬ 
tive-time  errors  at  any  station  pair  within  the  network. 
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All  of  the  time  data  used  in  this  study  were  derived 
from  nuclear  explosions  detonated  within  the  Nevada  Test  Site 
(NTS)  area  (Figure  1).  Table  I  lists  the  event  information. 

The  events  selected  as  references  from  which  residuals  and 
anomalies  were  measured  are  Bilby  or  Tan.  the  series  of  Bronze. 
Corduroy,  and  Buff,  or  the  series  of  Nash,  Agile,  and  Commodore. 
These  several  events  were  necessary  so  that  residuals  or  anoma¬ 
lies  could  be  obtained  for  an  adequate  number  of  recording 
stations.  Table  II  lists  the  station  anomalies,  relative  to 
RK-ON,  computed  from  the  travel-time  tables,  in  addition  to  the 
residuals  computed  from  the  Herri.n/1961  table.  A  key  letter, 
indicating  which  series  of  events  were  used  as  references  in 
determining  the  corrections  is  given  in  the  last  column. 

Station  records  of  all  explosions  were  routinely  read, 
with  the  identical  networks  of  stations  and  arrival  times  being 
used  for  both  programs.  Depths  for  all  events  were  restrained 
to  the  surface.  Raw  arrival  times  (Table  III)  were  input  to 
SHIFT  and  used  within  the  program  in  conjunction  with  the 
appropriate  input  travel- time  anomalies.  For  LOCATE,  station 
arrival  times  were  corrected  by  the  residuals  prior  to  their 
input  into  the  program.  These  input  times  for  LOCATE  are  ob¬ 
tained  by  subtracting  the  residuals  in  Table  II  from  the  arrival 

times  input  to  SHIFT  (Table  III)  . 

A  description  of  the  use  of  program  SHIFT  is  given  in 

Appendix  I  of  the  original  report. 

without  Time  Corrections  -  Table  IV  lists  the  location 

errors,  in  kilometers,  when  neither  residuals  nor  anomalies  are 
applied  to  the  event  arrival  times  (the  times  have  been  corrected, 
however,  for  station  elevation  and  ellipticity) .  The  mean  error 
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TABLE  I 


Event  Information 


Event 

Date 

Ori 

gin 

Time 

Latitude 

Longitude 

AUK 

02 

Oct 

64 

20 

03 

00.0 

37.078N 

116. 009W 

FORE 

16 

Jan 

64 

16 

00 

00.1 

.142 

.  049 

DUMONT 

19 

May 

66 

13 

56 

28.1 

.111 

.  058 

CHARTREUSE 

06 

May 

66 

15 

00 

00.1 

.  348 

.  322 

TURF 

24 

Apr 

64 

20 

10 

00.2 

.  150 

.  055 

KLICKITAT 

20 

Feb 

6  4 

15 

30 

00.1 

.151 

.  040 

PTLEDKIVER 

02 

Jun 

66 

15 

30 

00.1 

.  227 

.  055 

BRONZE 

23 

Jul 

65 

17 

00 

00.0 

.  098 

.  033 

CORDUROY 

03 

Dec 

65 

15 

13 

02.1 

.165 

.052 

BUFF 

16 

Dec 

65 

19 

15 

00.0 

.073 

.  029 

GREELEY 

20 

Dec 

66 

15 

30 

00.1 

.  302 

.  408 

PIRANHA 

13 

May 

66 

13 

30 

00.0 

.087 

.034 

NASH 

19 

Jan 

67 

16 

45 

00.1 

.  144 

.  136 

BOURBON 

20 

Jan 

67 

17 

40 

04 . 1 

.100 

.  004 

AGILE 

23 

Feb 

67 

18 

SO 

00.0 

.  127 

.  067 

COMMODORE 

20 

May 

67 

15 

00 

00.2 

.  130 

.064 

SCOTCH 

23 

May 

67 

14 

00 

00.0 

.  275 

.  370 

BILBY 

13 

Sep 

63 

17 

00 

00.1 

.061 

.  022 

TAN 

03 

Jan 

66 

14 

00 

00.0 

.  068 

.  035 

No .  of 
St a t ions 

6 

]  2 
8 
7 

10 

9 

7 
9 
9 

8 

13 
9 
7 

4 

7 

5 
5 

14 

8 


TABLE  II 

Station  Anomalies  and  Residuals 


Anomaly  Relative  to  RK-ON 


Residual 


Event 
* 


Station 

Herrin  61 

Herrin  66 

JB 

Herrin  61 

Key 

AD- IS 

+1.11  sec 

+0.04  sec 

+0.57 

sec  -1.1  sec 

B 

AX2AL 

+2.21 

+1.61 

+2.23 

+0.2 

A 

BE-FL 

+1.11 

+1.11 

+1.15 

-0.9 

A 

BL-WV 

+0.66 

+0.61 

+0.84 

-1.6 

A 

BR-PA 

+0.74 

+0.76 

+0.84 

-1.5 

A 

CP008 

+1.51 

+0.86 

+1.55 

-0.6 

A 

DH-NY 

+0.09 

+0.03 

+0.11 

-2.2 

A 

EB-MT 

-0.84 

+0.08 

-0.68 

-3.1 

A 

EN-MO 

+0.08 

+0.54 

+0.04 

-2.0 

B 

EU-AL 

+2.66 

+1.94 

+2.71 

+0.6 

A 

GG-GR 

+1.06 

+1.00 

+2.07 

-1.2 

A 

HN-ME 

+1.11 

+0.55 

+0.97 

-1.0 

A 

KC-MO 

+1.71 

+1.81 

+1.45 

-0.7 

D 

LV-LA 

+1.51 

+1.90 

+1.48 

-0.6 

A 

LZ-BV 

+0.90 

+0.94 

+1.19 

-1.9 

A 

NP-NT 

+1.88 

+0.96 

+1.52 

-0.2 

A 

OO-NW 

-0.01 

+0.01 

+0.52 

-2.4 

A 

PG-BC 

+2.66 

+2.94 

+2.57 

+0.3 

C 

PZ-PR 

+2.28 

+1.17 

+1.53 

+0.1 

A 

SI-BC 

+1.87 

+2.81 

+2.05 

-0.2 

A 

SV2QB 

+0.35 

-0.37 

+0.16 

-2.0 

A** 

SV3QB 

+0.35 

-0.37 

+0.16 

-2.0 

A 

WH2YK 

+1.36 

+1.13 

+1.50 

-0.8 

C 

RK-ON 

0 

0 

0 

-2.1 

A 

*A  =  Bilby-Tan 
B  =  Bronze-Corduroy- Buff 
C  =  Nash- Agile-Commodore 
D  =  Greeley 

**Set  equal  to  anomaly  for  SV3QB  determined  from  A 


TABLE  III  -  Arrival  Time  Data 


I 
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TABLE  IV 

Location  Errors  when  Neither  Residuals 
nor  Anomalies  are  Used 

NO.  OF 


•EVENT 

STATIONS 

L0CATE-H61 

SHIFT-H61 

SHIFT-H66 

SHIFT 

Buff 

8 

9.1 

8.6 

10.9 

5.7 

Turf 

10 

14.5 

15.0 

8.0 

18.3 

Corduroy 

9 

11.6 

10.8 

5.2 

6.3 

Nash 

7 

17.2 

18.1 

13.4 

23.1 

Bourbon 

4 

48.0 

49.0 

17.6 

42.6 

Piledriver 

7 

39.  6 

39.3 

37.7 

41.9 

Bronze 

9 

13.3 

13.1 

8.5 

12.3 

Auk 

6 

7.1 

7.2 

3.7 

4.8 

Piranha 

9 

41.1 

40.4 

45.5 

46.3 

Fore 

12 

38.1 

38.5 

20.0 

36.6 

Greeley 

13 

10.1 

12.3 

11.6 

25.5 

Dumont 

8 

43.0 

43.6 

44.9 

47.1 

Chartreuse 

7 

60.0 

59.8 

48.5 

62.8 

Commodore 

5 

11.1 

10.9 

23.3 

12.6 

Agile 

7 

23.4 

23.0 

15.0 

27.7 

Scotch 

5 

17.4 

17.8 

26.6 

19.6 

Klickitat 

9 

33.5 

33.7 

12.8 

38.2 

I 

438.1 

441.1 

353.2 

471.4 

Mean  error,  km 

25.8 

25.9 

20.8 

27.7 

for  17  events  using  LOCATE  is  seen  to  be  25.8  km  and  using  SHIFT 
(Herrin  61)  25.9  km/  the  results  from  the  two  programs  using  the 
same  travel-time  table  are  essentially  in  agreement.  With  the 
J-B  and  Herrin  66  tables,  the  results  using  SHIFT  are  27.7  km 
and  20.8  km  respectively.  It  is  not  known  at  this  time  if  the 
apparently  better  results  obtained  by  using  the  Herrin  66  table 
are  significant  or  not.  It  would  be  necessary  to  locate  a 
larger  sample  of  events  to  determine  the  effectiveness  of  this 
table  for  the  networks  used  in  this  study  and  for  the  NTS  area. 

The  resultant  average  errors  of  about  26  km  demonstrate 
the  best  one  could  hope  to  do  when  no  allowances  are  made  for 
residuals  (LOCATE)  or  anomalies  (SHIFT). 

With  Time  Corrections  —  When  the  residuals  or  anomalies 
are  determined  for  a  particular  region,  such  as  the  Nevada  Test 
Site,  and  for  each  station  which  is  to  be  used  in  subsequent 
location  networks,  the  resultant  location  errors  can  be  reduced 
by  at  least  an  order  of  magnitude.  Table  V  lists  the  location 
errors,  in  kilometers,  when  either  residuals  (LOCATE)  or  rela- 
tive  anomalies  (SHIFT)  are  applied  to  the  event  arrival  times. 

The  mean  error  for  17  events  using  LOCATE  is  now  seen  to  be 
2.98  km  and  using  SHIFT  (Herrin  61)  2.86  km.  Again,  the  mean 
values  from  either  program  are  in  agreement,  although  individual 
event  locations  differ  by  as  much  as  4.1  km  between  the  two 
methods.  This  difference  implies,  of  course,  that  the  programs 
are  computing  in  significantly  different  ways  but  yield  about 
the  same  answer  on  the  average.  Using  the  J-B  and  Herrin  66 
tables,  the  SHIFT  mean  errors  are  2.92  km  and  2.59  km,  respect¬ 
ively,  which  suggests  that  when  anomalies  are  applied,  any 
reasonable  travel-time  table  is  adequate. 
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TABLE  V 

Location  Errors  When  Residuals 
or  Anomalies  are  Used 

NO.  OF 


EVENT 

STATIONS 

L0CATE-H61 

SHIFT-H61 

SHIFT-H66 

SHIFT- 

Buff 

8 

2.4 

0.1 

0.3 

0.2 

Turf 

10 

0.  5 

0.1 

0.2 

0.1 

Corduroy 

9 

1.6 

1.0 

0.9 

0.9 

Nash 

7 

5.8 

1.7 

1.4 

1.8 

Bourbcn 

4 

2.8 

1.9 

2.1 

1.9 

Piledriver 

7 

2.0 

2.9 

1.5 

3.3 

Bronze 

9 

4.2 

3.0 

2.7 

2.9 

Auk 

6 

3.1 

3.0 

2.9 

3.0 

Piranha 

9 

0.8 

3.1 

2.7 

3.2 

Fore 

12 

3.3 

3.2 

2.9 

3.1 

Greeley 

13 

1.0 

3.3 

3.8 

3.3 

Dumont 

8 

2.2 

3.4 

3.0 

3.4 

Chartreuse 

7 

5.7 

3.6 

3.7 

4.2 

Commodore 

5 

4.1 

3.7 

3.2 

3.3 

Agile 

7 

1.5 

3.9 

2.4 

3.1 

Scotch 

5 

3.S 

5.0 

4.8 

5.7 

Klickitat 

9 

6.1 

6.0 

5.6 

6.2 

l 

50.7 

48.9 

44.1 

49.6 

Mean  error, 

km 

2.98 

2.86 

2.59 

2.92 

Confidence  regions  were  computed  within  SHIFT  in  the 
usual  manner  (e.g.,  Flinn,  1965)  for  all  events,  both  with  and 
without  travel-time  anomalies.  Table  VI  lists,  by  event,  the 
areas  of  the  computed  ellipses  and  the  factors  of  reduction  in 
ellipse  areas  when  anomalies  are  used.  Without  anomalies,  two 
events  (Fore  and  Chartreuse)  were  not  within  the  ellipse  at  the 
95%  level.  With  anomalies,  all  events  were  within  the  ellipse 
and  the  ellipses  were  reduced  in  area  by  factors  of  5  to  152, 
with  an  average  reduction  of  about  45.  This  reduction  in  el¬ 
lipse  area  points  out  the  necessity  and  value  of  travel-time 
anomalies. 

The  dashed-line  polygon  in  Figure  2  is  the  contour  of 
the  (known)  maximum  relative  error  (0.18  sec)  at  the  true  epi¬ 
center  marked  at  X.  The  circle  enclosing  the  value  0.11  is  the 
maximum  relative  error  as  a  result  of  the  final  solution  ob¬ 
tained  with  SHIFT-H61.  The  solid-line  polygon  is  the  contour 
of  the  maximum  relative  error  estima ted  to  be  0.4  sec  (0.3  sec 
higher  than  that  at  the  final  solution)  for  this  event.  A 
seismic  analyst  can  usually  estimate  his  reading  errors  quite 
well,  and  if  the  effects  of  previously  determined  relative 
travel-time  anomalies  are  removed  (they  must  be  determined, 
not  estimated.)  ; the  estimates  can  be  used  to  contour  the  maximum 
relative  error.  A  relative  error  estimate  of  0.4  sec  is  liberal 
for  an  event  of  the  size  of  Bourbon,  but  it  is  an  example  of 

the  manner  of  using  any  estimate.  The  approximate  area  of  the 

2 

estimated  polygon  is  less  than  340  km  .  The  area  of  the  cor- 

2 

responding  ellipse  previously  computed  (Table  VI)  is  7118  km  . 

Subroutine  SI^RID  produces  an  output  of  network  stan¬ 
dard  deviations,  in  addition  to  the  maximum  relative  errors, 
for  a  similar  set  of  grid  positions.  Figure  3  shows  the  zero- 
mean  standard  deviation  output  for  Bourbon,  with  anomalies. 
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TABLE  VI 


Confidence  Regions 
SHIFT-61 


Event 

Ellips 
Without  anomalies 
a 

e  Areas  kni 

With  anomalies 
b 

Improvement 
Factor,  a/b 

No .  of 
Stations 

Fore 

2  0  5  9  * 

348 

6 

12 

Buff 

2402 

50 

48 

8 

Chartreuse 

1041* 

131 

8 

7 

Auk 

10035 

66 

152 

6 

Bourbon 

206088 

7118 

29 

4 

Dumont 

1495 

52 

29 

8 

Agile 

4020 

234 

17 

7 

Nash 

7256 

113 

64 

7 

Commodore 

27698 

288 

96 

5 

Greeley 

2620 

96 

27 

13 

Klickitat 

2384 

202 

12 

9 

Turf 

2513 

524 

¥  5 

10 

Piranha 

2159 

15 

144 

9 

Scotch 

37869 

2309 

16 

5 

Corduroy 

1564 

25 

62 

9 

Bronze 

2038 

25 

81 

9 

Piledriver 

1593 

79 

20 

7 

*Conf idence 

ellipse  does  not 

contain  true  epicenter. 
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The  following  conclusions  are  made  concerning  the  results 
obtained  by  analyzing  19  nuclear  explosions  detonated  within  a 
2500  km2  area  of  the  Nevada  Test  Site.  Seismograms  were  used 
from  teleseismic  stations  forming  networks  of  four  to  thrrteen 

stations. 

1.  For  limited-station  teleseismic  networks,  the  loca¬ 
tion  capability  for  17  events  occurring  in  the  Nevada  Test  Site 
and  without  applying  previously-determined  residuals  or  travel¬ 
time  anomalies,  is  about  26  km,  regardless  of  the  program  used 
(LOCATE  or  SHIFT)  and  regardless  of  the  travel-time  table  em¬ 
ployed  (Herrin  61,  Herrin  66,  or  J-B  tables). 

2.  The  location  capability  for  the  same  networks  and 

time  data  is  better  than  3  km  when  previously-determined  resid¬ 
uals  or  travel-time  anomalies  are  applied,  regardless  of  program 
and  travel-time  table.  That  is  to  say,  when  these  travel-time 
curves  can  equivalently  be  replaced  by  observed  travel  times 
from  each  station  to  an  accurately  known  explosion  point,  the 
particular  curve  selected  as  a  standard  is  essentially  irrele¬ 
vant  and  epicenters  of  other  nearby  explosions  can  be  located 

within  the  indicated  error  limit. 

3.  The  effect  of  the  number  of  stations  in  the  range 

from  3  or  4  to  13  on  the  location  capability  of  the  networks 
used  in  the  study  is  negligible,  either  with  or  without  travel- 

time  anomalies. 

4.  The  effect  of  the  range  of  epicentral  distances 
(distance  aperture)  on  the  capability  of  locating  known  sur¬ 
face  events  with  the  same  networks  is  negligible,  either  with 
or  without  travei-time  anomalies. 
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5.  The  azimuth  aperture  has  an  observable  effect  on  the 
location  capability  of  the  networks  where,  for  apertures  down  to 
60°,  the  location  errors  without  travel-time  anomalies  are  as 
high  as  60  km,  and  with  anomalies  6  km. 

6.  The  areas  of  computed  confidence  ellipses  can  be 
reduced  by  factors  of  1/5  to  1/152  with  the  application  of  travel¬ 
time  anomalies  and  still  enclose  the  event  locations  which,  for 
this  study,  are  accurately  known,  as  are  the  explosions  used  to 
calibrate  the  particular  epicenter-station  paths  used  in  the 
study.  In  effect,  these  reduced  ellipses  represent  uncertain¬ 
ties  in  the  position  of  the  epicenters  relative  to  the  position 

of  the  calibration  events. 

7.  With  as  few  as  4  or  5  stations,  computed  confidence 

ellipses  with  anomalies  are  unrealistically  large  due  to  the 

small  number  of  statistical  degrees  of  freedom.  (The  Bourbon 

2 

explosion,  with  4  stations,  has  an  ellipse  7,000  km  and  a  loca¬ 
tion  error  of  1.9  km).  Estimates  of  the  maximum  relative  time 
errors  for  a  network  of  stations  permits  more  realistic  confi¬ 
dence  limits  to  be  set  on  the  solution.  In  this  way.  Bourbon's 

2 

confidence  polygon  is  about  340  km  .  Also,  as  few  as  three 
stations  can  be  used  to  obtain  the  confidence  limits  when  maxi¬ 
mum  relative  errors  are  estimated. 

8.  Travel-time  anomalies,  computed  from  a  few  selected 

events,  show  some  variability  across  the  NTS  area  in  Figure  1, 

but  the  effect  on  relative  location  accuracy  is  small.  The  size 

2 

of  the  area  in  Figure  1  is  about  2500  km  ,  implying  that  fairly 
large  regions  are  involved  for  determining  reasonably  constant 
anomalies,  so  the  problem  of  time  calibrating  a  stationary  net¬ 
work  for  any  region  of  the  earth  from  which  either  explosions 
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at  accurately  known  locations  or  earthquakes  (bias  effects  not 
included)  are  recorded  is  not  formidable.  In  either  case,  if 
accurate  locations  are  not  independently  known,  the  epicentral 
solutions  with  anomalies  included  reduce  to  locations  relative 
to  one  another  with  the  actual  error  of  the  whole  set  remaining 
unknown. 

9.  Relative  anomalies  for  earthquake  regions  can  be 
simpler  to  assess  than  absolute  residuals  because  (a)  origin 
time  errors  are  eliminated;  (b)  first-extremum  anomalies  and 
first  motion  anomalies  from  several  events  can  be  combined 
(except  when  obvious  period  differences  are  noted) ;  and  (c) 
depending  on  the  region  and  on  the  network  geometry.-  the  regions 
for  which  the  calibrations  need  to  be  determined  can  be  fewer 
in  number  and  larger  in  area. 
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B.  Preliminary  Beamforming  Study  of  the  TFO-37  Array 

This  study  describes  a  preliminary  evaluation  of  the 
new  TFO  array  in  terms  of  the  amount  of  signal  loss,  noise 
attenuation,  and  S/N  gain  which  result  from  beamforming  digital 
recordings  of  four  teleseismic  events.  The  array  consists  of 
37  JM  short-period  vertical  seismometers  arranged  to  form  three 
concentric  hexagons  containing  six,  twelve,  and  eighteen  sensors, 
respectively,  with  an  additional  sensor  in  the  center.  The  ap¬ 
proximate  diameter  of  the  outside  hexagon  is  thirty  kilometers, 
and  intersensor  spacing  is  approximately  five  kilometers. 

In  addition  to  presenting  the  results  of  a  basic  in¬ 
vestigation,  another  intention  is  to  assemble  information  per¬ 
taining  to  the  array  in  one  report  which  can  serve  as  a  reference 


for  future  studies.  With  this  in  mind,  the  following  appendices 

are  included  in  the  original  report: 

Appendix  A:  Map  of  Array 

Seismometer  description  and  locations 
in  latitude  and  longitude 

Digital  tape  format 

Digital  programs  written  for  the  evaluation 

Calibration  parameters 

Sample  output  from  which  Table  VII  was 
prepared 

JM  frequency  response 
The  LASA  Bulletin  statistics  for  the  data  are  shown  in 
Table  VIII.  Since  this  report  includes  only  a  small  number  of 
events,  one  from  each  quadrant  was  selected  to  given  full  azimuthal 
coverage^ 


Appendix  B: 

Appendix  C: 
Appendix  D: 
Appendix  E: 
Appendix  F: 

Appendix  G: 
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Table  VII  summarizes  the  final  results  of  the  amplitude 
analysis  and  contains  three  major  headings.  Signal,  Noise,  and 
S/N.  As  shown  in  Table  VII,  there  is  a  slight  loss  in  signal 
amplitude  between  the  mean  and  the  sum  possibly  due  to  small 
misalignments.  In  addition,  the  apparent  db  loss  also  can  be 
attributed  to  the  signal/noise  ratio.  The  Tunisia  earthquake 
is  poorly  recorded  and  the  “P"  probably  contains  some  noise 
which  precludes  an  accurate  estimate  of  signal  amplitude  during 
the  phasing  process.  A  deterioration  in  amplitude  (my  )  is  ob¬ 
vious  between  the  SDL  and  Fc  filter  (Figure  4)  amounting  to 
about  4  db  in  all  cases.  This  suggests  that  the  SDL  filter  is 
more  effective  in  preserving  signal  amplitude  within  the  fre¬ 
quency  range  of  the  events  studied.  On  the  other  hand,  an  exam¬ 
ination  of  Figure  5  reveals  the  undesirable  loss  of  first  motion 
as  well  as  a  loss  of  signal  character  caused  by  the  SDL  filter 
which  appears  to  over-simplify  the  signature. 

With  respect  to  the  noise,  both  prefilters  yield  an 
acceptable  reduction  in  rms,  but  the  Fc  is  the  better  of  the  two. 
This  is  due  to  the  removal  of  the  low  frequency  microseisms 
which  predominate.  The  amplitude  of  the  maximum  noise  excursion, 
indicated  by  the  heading  RANGE,  is  also  shown.  The  ratio  of 
the  Range/rms  is  approximately  3/1,  a  value  normally  considered 
when  referring  to  peak  versus  average  rms  noise. 

The  last  heading  in  the  table  concerns  S/N  ratio. 

Again,  both  filters  function  well  but  the  Fc  responds  best  by 
about  3  db  over  the  SDL. 
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Figure  4. 


Bandpass  filters  used  in  data  preparation. 


Power  spectra  were  completed  for  the  100  seconds  of 
noise  preceding  each  of  the  four  events  by  a  digital  program 
based  on  the  Blackman  and  Tukey  method.  Figures  7  through  10 
each  illustrate  the  results  of  spectral  analysis  of  noise 
occurring  prior  to  the  signal  on  the  phased  sum  and  the  average 
noise  spectra  of  outputs  recorded  by  18  to  19  individual  elements. 
It  was  felt  that  taking  the  average  of  a  number  of  sensors  was 
better  than  selecting  just  one,  and  that  spatial  sampling  of  the 
array  would  give  about  the  same  results  as  proces&i.ng  the  total. 

The  spectral  shape  is  approximately  the  mime  for  any 
given  filter  and  will  not  vary  much  whether  the  analysis  is  per¬ 
formed  on  an  individual  trace  or  the  phased  sum  of  many.  In  all 
examples  shown,  summing  the  elements  attenuates  the  noise  at, 
for  example,  1  cps  approximately  15  db  over  the  average  single 
instrument  which  is  equal  to  a  factor  of  N.  The  best  S/N  ratio, 
again  referring  to  Table  VII,  is  several  db  lower.  For  com¬ 
parative  purposes.  Figure  6  shows  the  results  of  raw  data  proces¬ 
sing  of  the  Hokkaido  event.  This  reveals  that  the  lower  fre¬ 
quencies  are  basically  coherent  across  the  array  and  are  pre¬ 
served  even  after  time  shifting. 

The  results  of  prefiltering  and  beamforming  TFO  digital 
recordings  of  four  teleseismic  events  indicate  that: 

1.  The  P  phase  from  strongly  recorded  events  is  atten¬ 
uated  only  slightly  by  beamforming.  However,  it  appears  that 
the  weaker  the  event,  the  greater  the  loss. 

2.  The  rms  of  the  noise  is  reduced  by  11-14  db, 
depending  upon  the  filter,  after  summing. 

3.  The  zero-to-peak  maximum  noise  excursion  in  a  100- 
second  window  is  about  three  times  that  of  the  rms  obtained 
from  the  same  sample. 
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Power  Soectra  of  unfiltered  noise 
Hokkai  J.O,  25  February  1968. 
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Figure  7.  Noise  power  spectra,  Hokkaido, 
25  February  1968. 
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Figure  9.  Noise  power  spectra,  Peru, 

24  February  1968. 
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Figure  10.  Noise  power  spectra,  Tunisia, 
25  February  1968. 


4.  There  is  also  an  11-14  db  improvement  in  S/N 
ratio  by  beamforming. 

5.  Beamforming  produces  about  15  db  noise  reduction 
at  1  cps  as  shown  by  the  power  spectra. 

6.  Neither  filter  used  is  completely  satisfactory. 

The  SDL  seems  to  preserve  signal  amplitude  at  the  expense  of 
character,  but  does  not  reduce  the  noise  to  a  minimum.  The  Fc 
filter  reduces  the  noise  almost  to  the  system  background,  but 
attenuates  the  signal  amplitude  while  preserving  the  first 
motion  and  shape. 

As  a  result  of  this  study,  three  areas  of  interest  are 
indicated  which  should  be  investigated;  they  are: 

1.  The  use  of  average  observed  travel-time  anomalies 
and  cross-correlations  to  reduce  beamforming  signal  loss. 

2.  The  design  and  evaluation  of  a  set  of  bandpass 
filters  to  determine  the  most  efficient  with  respect  to  signal 
preservation  and  noise  reduction. 

3.  The  comparison  of  TFO-37  with  LASA,  NORSAR  and  the 
other  arrays  in  terms  of  their  respective  detection  and  identi¬ 
fication  capabilities. 

c*  P  Wave  Parameters  Measured  at  the  Montana  LASA 

The  observed  times  of  first  arrival  P-waves  recorded 
by  a  horizontal  array  of  seismometers  such  as  the  one  at  LASA 
provide  measurements  that  can  be  used  to  compute  (1)  the  recip¬ 
rocal  horizontal  phase  velocity  (p) ,  and  (2)  the  deviation  of 
direction  of  arrival.  The  raypath  parameter  p  is  the  magnitude 
of  the  slope  of  the  travel— time  curve,  and  the  azimuth  deviation 
is  the  angle  between  the  back  azimuth  to  the  epicenter  and  the 
apparent  arrival  direction.  Both  of  these  wavenumber  parameters 
yield  information  about  the  crust  and  upper  mantle  structure 
beneath  the  array  and  the  variation  of  velocity  in  the  mantle. 
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The  first  step  in  computing  these  quantities  is  to  fit 
a  plane  wave  front  to  the  first  arrival  P-waves.  This  was  done 
by  the  least-squares  method  described  in  the  original  report. 
Making  the  assumption  that  the  wave  front  is  planar  imposes 
the  following  practical  limitations  on  events  studied:  (1) 
the  events  which  can  be  considered  must  occur  at  distances 
greater  than  2900  km  (26°)  from  the  array;  (2)  the  maximum 
dimension  of  the  array  must  be  less  than  100  km  as  measured 
from  the  reference  station. 

The  errors  due  to  scatter  in  observed  times  can  be 
reduced  to  a  minimum  by  using  an  array  with  the  largest  pos¬ 
sible  diameter.  The  F-ring  of  subarrays  at  LASA  has  a  diameter 
of  200  km,  which  is  ideally  suited  for  this  purpose. 

The  raypath  parameter,  p,  and  the  azimuth  deviation 
have  been  computed  for  over  600  earthquakes  recorded  at  LASA 
during  its  first  two  years  of  operation.  The  epicenters  of 
T47  events  are  located  along  a  profile  extending  to  the  north¬ 
west  in  the  range  of  azimuth  300°-320°.  These  are  so  distrib¬ 
uted  as  to  provide  essentially  continuous  coverage  of  the 

o  o 

distance  range  26  -97  . 

To  reduce  scatter  in  p  and  the  azimuth  deviation  we 
consider  only  171  events  which  satisfy  the  following  criteria: 

(1>  the  event  must  have  been  clearly  recorded  at  the  reference 
station,  A0,  and  at  all  four  of  the  stations  in  the  F-ring  of 
subarrays  (because  the  determination  of  the  least-squares 
source  of  an  event  located  within  this  range  of  azimuth  criti¬ 
cally  depends  upon  the  delay  times  at  subarrays  F2  and  F4) ;  (2) 
the  event  must  have  been  recorded  at  least  at  eight  of  the  other 
16  subarrays. 
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Figure  11  shows  the  variation  with  distance  of 
measured  azimuth  deviations  along  the  northwest  profile.  The 
two  important  features  of  the  observed  variation  are:  (1) 
the  average  value  of  the  azimuth  deviations  for  all  events 
along  this  profile  is  -0.5°,  and  (2)  the  variation  in  the 
azimuth  deviations  reflects  the  change  in  epicentral  azimuth 
of  these  events.  This  becomes  apparent  when  we  examine  the 
azimuth  deviations  as  a  function  of  azimuth  along  the  profile. 

For  the  northwest  profile  the  opicents  in  the  distance 

range  32°  and  47°  are  within  the  azimuth  range  300°  -  305°.  All 

but  four  of  the  corresponding  azimuth  deviations  have  values 

between  -0.5°  and  -2.0°.  In  the  distance  range  52°  -  75  all 

o  o 

but  one  of  the  epicents  are  between  the  azimuths  310  and  317  , 
and  all  but  three  of  the  azimuth  deviations  corresponding  to 
these  events  are  between  -0.5°  and  +1.1°.  In  addition,  the 
value  of  the  azimuth  deviations  changes  most  in  the  distance 
intervals  26°-32°,  47°-52°,  and  75°-97°  where  the  corresponding 
epicentral  azimuths  are  undergoing  the  most  change. 

Figure  12  shows  the  variation  of  p  with  distance  for 
events  in  the  northwest  profile.  The  observations  provide 
almost  complete  coverage  between  26°  and  97°.  With  the  ex¬ 
ception  of  the  abrupt  discontinuity  in  the  curve  at  87  ,  all 
but  five  of  the  observed  values  of  p  are  with  +0*05  seconds 
per  degree  of  the  smooth  curve  fitted  to  the  average  value. 

The  discontinuity  which  occurs  at  87°  is  discussed  below. 

This  curve  can  be  compar  ed  with  the  similar  curve 
which  Chinnery  and  Toksdz  (1967)  constructed  from  a  smaller 
number  of  events  along  the  same  profile.  The  curve  that  we 
present  includes  171  events  for  which  the  delay  times  at  all 
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four  stations  in  the  F-ring  of  subarrays  are  available;  Chinnery 
and  Toksttz  (1967,  p.  204)  included  only  78  events  on  their  curve. 
Our  data  are  more  complete  in  the  intervals  between  32  -40  , 
5]L0_550f  and  72°-78°,  which  are  of  particul  ar  interest. 

Figure  13  shows  the  comparison  between  a  smooth  curve 
fitted  through  our  observed  points  and  the  corrected  observa¬ 
tions  given  by  Chinnery  and  Toksttz  (1967,  Table  1,  p.  215-216). 

For  reference.  Figure  13  also  shows  the  Jef freys-Bullen  curve. 

Our  curve  includes  only  small  adjustments  required  for  smoothing. 
The  Chinnery-Toksttz  curve  includes  a  constant  correction  of  +0.05 
seconds  per  degree  to  make  the  observations  consistent  with  the 
absolute  travel-times  (see  Chinnery  and  Toksttz,  1967,  p.  214). 

The  Jef freys-Bullen  curve  was  computed  by  fitting  a  7-point 
polynomial  to  the  J-B  travel-times  (see  Chinnery  and  Toksttz,  1967, 

p.  204) . 

Comparison  of  the  two  observed  curves  shows  three  intervals 
in  which  Chinnery  and  Toksttz  (1967)  extrapolated  their  values 
(32°-40°,  51°-55°,  and  72°-78°) .  The  additional  points  on  our 
curve  demonstrate  that  the  curve  actually  lies  below  the  position 
that  they  suggest.  In  the  intervals  32°-40°  and  72  -78  ,  where 
their  curve  is  projected  above  the  corresponding  portions  of  the 
j-B  curve,  the  new  data  establish  that  the  observed  curve  follows 
these  portions  of  the  J-B  curve  much  more  closely.  For  the  in¬ 
terval  between  52°-55°,  where  Chinnery  and  Toksttz  place  their 
curve  only  slightly  below  the  J-B  curve,  our  additional  points 
suggest  that  this  portion  of  the  observed  curve  falls  well  be¬ 
low  the  Jef freys-Bullen  curve.  The  remaining  portions  of  our 
curve,  out  to  a  distance  of  87°,  are  in  good  agreement  with 
Chinnery  and  Toksttz,  except  for  the  correction  of  +0.05  seconds 
per  degree.  This  correction  has  not  been  added  to  our  curve. 
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The  segments  along  our  curve  which  differ  from  the 
Chinnery-Toksttz  curve  are  of  particular  interest.  Chinnery 
and  Toksdz  (1967,  p.  233)  conclude  from  their  results  that 
the  points  where  the  slope  increases  abruptly  (35  ,  52  and 
70°  on  their  curve)  reflect  points  of  inflection  in  the 
velocity-depth  function.  They  further  conclude  that  these 
points  of  inflection  occur  at  depths  of  approximately  tOO,  1300, 
and  2000  kilometers,  respectively.  The  new  detail  provided  by 
our  values  of  p  shows  abrupt  slope  increases  at  32  ,  51  ,  58  , 
and  71°.  Our  data  thus  indicate  that:  (1)  on  the  northwest 
profile  an  additional  point  of  inflection  occurs  in  the  velocity- 
depth  function  at  a  depth  of  approximately  1500  kilometers,  and 
(2)  the  upper  point  of  inflection  occurs  at  a  depth  of  approxi¬ 
mately  750  kilometers. 

Because  our  values  of  p  were  determined  using  only 
first  arrivals  they  do  not  show  any  travel-time  triplications 
directly,  but  they  do  show  two  features  expected  at  triplication 
distances.  These  are:  (1)  relative  absence  of  points  where  the 

slope  increases  rapidly,  which  may  result  from  the  reduced  am¬ 
plitude  of  the  P  wave  arrival  where  triplication  occurs;  (2) 
relative  abundance  of  points  in  the  interval  immediately  pre¬ 
ceding  the  c'  ange  in  slope,  as  a  result  of  focusing  by  the 
bottom  of  the  velocity  transition  layer.  Our  curve  shows  a 
concentration  of  points  in  the  intervals  29°-31°,  49  -51  ,  56  - 
57°,  and  68°-71°  and  a  corresponding  lack  of  points  in  the  in¬ 
tervals  31°-33°f  51°-53° ,  57°-58° ,  and  71°-72°. 

The  abrupt  discontinuity  in  the  curve  of  p  vs.  distance 
for  the  northwest  profile  at  87°  is  probably  due  to  a  large 
change  in  epicentral  azimuth  of  the  corresponding  events.  All 
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of  the  p  values  in  the  interval  83°-88°  fall  on  the  lower  port.io 
of  the  curvo  and  the  corresponding  epicenters  are  in  the  range 
of  azimuth  300°-308°.  The  epicentral  azimuths  of  all  the  events 
in  the  distance  interval  87°-97°,  where  the  p  values  fall  on  the 
upper  portion  of  the  curve,  are  between  309°  and  314°.  Thus, 
the  systematic  change  in  the  p-distanc<?  curve  reflects  a  syste¬ 
matic  change  in  the  corresponding  epiceniral  azxmuths  along  the 
northwest  profile. 

The  southeast  profile  includes  122  events  in  the  range 
of  azimuth  140°-160°.  forty  other  events  whose  epicenters  lie 
on  the  profile  were  omitted  because  one  or  more  stations  in  the 
F-ring  did  not  record  the  signal.  The  122  events  finally  in¬ 
cluded  in  the  southeast  profile  satisfy  all  of  the  criteria  for 
quality  stated  in  the  previous  section. 

Figure  14  shows  the  variation  with  distance  of  measured 
azimuth  deviations  along  the  southeast  profile.  Two  important 
observations  about  azimuth  deviations  along  this  profile  are: 

(1)  the  average  value  of  +1.0°,  and  (2)  a  strong  relationship 
with  epicentral  azimuths. 

In  the  distance  interval  38°-72°  all  of  the  epicenters 
are  in  the  range  of  azimuth  142°-150°  and  all  but  three  of  the 
related  values  for  azimuth  deviations  are  between  +0.9°  and 
+2.5°.  In  the  distance  intervals  31°-38°  and  72°-97°,  where 
the  epicentral  azimuths  change  by  more  than  10° ,  the  corres¬ 
ponding  azimuth  deviations  also  undergo  the  greatest  rate  of 
change. 

Comparision  of  Figure  14  with  Figure  11  shows  that  the 
azimuth  deviation  along  the  southeast  profile  is  greater  than 
along  the  northwest  profile  by  at  least  1.5°.  This  demonstrates 
a  strong  function  of  azimuth  along  both  profiles. 
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Figure  14.  Observed  azimuth  deviations  for  events  with 

epicenters  located  along  the  southeast  profile, 
140°-160° . 


Figure  15  shows  the  variation  of  p  with  distance  for 
events  along  the  southeast  profile.  The  observed  points  provide 
good  coverage  of  the  distance  range  between  31°  and  86°  with 
the  exception  of  the  intervals  between  39°-42°,  43°-45°,  47°-52°, 
56°-60°,  and  63°-66°.  Five  observed  points  in  the  interval 
86°-97°  establish  the  approximate  position  of  this  portion  of  the 
curve.  All  but  four  of  the  observed  points  along  the  southeast 
profile  are  within  j+0  ,08  seconds  per  degree  of  the  smooth  curve 
fitted  to  the  average  value. 

The  lack  of  observed  data  in  the  above  intervals  creates 
some  doubt  as  to  the  actual  position  of  the  curve  in  these  inter¬ 
vals.  However,  available  data  do  show  the  following  features 
along  the  southeast  profile:  (1)  An  abrupt  increase  in  slope  or 
an  actual  triplication  is  clearly  indicated  in  the  distance  in¬ 
terval  34°-35°;  (2)  Similar  abrupt  increases  in  slope  or  tripli¬ 
cations  are  suggested  in  the  distance  intervals  49°-54°  and  63°- 
66°;  (3)  The  distance  interval  39°-42°  may  include  such  an 
abrupt  increase  in  slope;  (4)  There  are  no  clear  indications  of 
abrupt  increases  in  slope  at  distances  32°,  58  ,  or  71°,  where 
increases  in  slope  were  observed  on  the  northwest  profile. 

Figure  16  shows  a  smooth  curve  fitted  to  p  observed 
along  the  southeast  profile,  together  with  the  observed  curve 
for  the  northwest  profile.  The  extrapolated  distance  intervals 
are  shown  by  broken  lines.  Comparison  of  the  two  curves  shows 
that  the  southeast  profile  curve  is  displaced  from  the  north¬ 
west  profile  curve  by  an  amount  0.48  seconds  per  degree,  on  the 
average.  The  position  of  the  curve  from  the  southeast  profile 
after  a  constant  correction  of  -0.48  seconds  per  degree  has 
been  applied  is  also  shown  in  Figure  16. 
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Comparison  of  the  northwest  profile  with  the  south¬ 
east  shows:  (1)  the  abrupt  increases  in  slope  which  are 
clearly  indicated  at  32  and  71  on  the  northwest  profile  are 
not  at  all  suggested  by  the  data  available  from  the  southeast 
profile;  (2)  there  are  abrupt  slope  increases  along  the  south¬ 
east  profile  in  the  intervals  34°-35°  and  63°-66°  which  are 
clearly  not  indicated  for  the  same  intervals  on  the  northwest 
profile;  (3)  it  is  possible  that  the  abrupt  increases  in 
slope  which  occur  at  51°  and  58  on  the  northwest  profile  also 
occur  at  the  same  distances  along  the  southeast  profile,  al¬ 
though  the  sparse  data  from  the  southeast  profile  suggest  that 

such  a  change  does  occur  between  53°  and  54  but  does  not 
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occur  between  56  and  62  . 

As  for  the  northwest  profile,  each  place  where  the  p 
versus  distance  curve  has  an  abrupt  increase  in  slope  reflects 
a  corresponding  increase  in  the  velocity-depth  function. 
Present  data  from  the  southeast  profile  do  not  provide  enough 
detail  in  the  critical  distance  intervals  to  determine  the 
actual  depths  at  which  these  velocity  transition  zones  occur. 
However,  these  data  strongly  suggest  that  three  such  velocity 
transitions  occur  beneath  the  southeast  profile,  at  corres¬ 
ponding  depths  of  approximately  800,  1350  and  1800  kilometers 
(these  depth  estimates  are  crude  interpolations  made  using 
the  figures  of  Chinnery  and  Tokstiz,  1967).  Similar  transition 
zones  have  been  previously  observed  in  the  upper  mantle  by 
many  other  workers  (for  example,  Archambeau  et  al,  1967; 
Johnson,  1967;  Anderson,  1967). 
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The  average  value  of  the  observed  azimuth  deviations 
along  the  southeast  profile  is  1.5°  greater  than  that  for  the 
northwest  profile,  and  the  observed  p  values  for  the  southeast 
profile  are  greater  than  those  for  the  northwest  profile  by 
nearly  0.5  seconds  per  degree.  In  addition,  the  observed  dis¬ 
continuity  at  87  in  the  p  curve  for  the  northwest  profile  is 
apparently  related  to  a  systematic  change  in  the  epicentral 
azimuths.  These  observations  show  that  both  wave-number 
parameters  change  considerably  with  changes  in  epicentral 
azimuth. 

Azimuthal  changes  of  this  kind  can  result  from  either 
dipping  layers  or  horizontal  velocity  gradients  within  the 
crust  and  upper  mantle  beneath  the  array.  From  a  study  of  the 
relative  station  corrections  at  LASA,  Fairborn  (1966)  has  con¬ 
cluded  that  dipping  layers  cannot  account  for  the  magnitude  of 
the  required  corrections  but  that  horizontal  velocity  gradients 
may  account  for  the  required  magnitude  if  they  extend  into  the 
upper  mantle,  ^hus  we  conclude  that  the  azimuthal  variation  in 
p  and  azimuth  deviation  is  probably  caused  by  horizontal  velocity 
gradients  beneath  LASA  which  extend  into  the  upper  mantle. 

The  observed  variation  with  distance  of  p  along  both 
the  northwest  and  southeast  profiles  shows  abrupt  increases 
in  the  slope  of  the  curve  at  several  places.  However,  the 
distances  at  which  these  changes  occur  are  not  the  same  along 
the  two  profiles.  The  additional  detail  provided  by  the  new 
data  from  this  study  for  the  northwest  profile  suggests  that 
abrupt  increases  in  the  slope  of  the  p  distance  curves  reflect 
several  previously  observed  velocity  increases  within  the  mantle. 
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For  the  northwest  profile,  our  data  indicate  that  four  such 
velocity  transition  zones  occur,  at  depths  of  approximately 
750,  1300,  1500,  and  2000  kilometers.  For  the  southeast 
profile  the  data  are  not  as  detailed,  but  they  suggest  that 
at  least  three  such  velocity  zones  occur  at  depths  of  approxi¬ 
mately  800,  1350,  and  1800  kilometers.  Thus,  the  new  data 
presented  here  demonstrate  that  the  velocity-depth  function 
for  the  region  extending  northwest  from  LASA  is  distinctly 
different  from  that  for  the  region  to  the  southeast  for  depths 

between  700  and  2000  kilometers. 

Departures  of  times  from  the  least- squares  plane  wave 

front  are  the  stations  deviations.  The  results  suggest 
additional  work  to  determine  if  station  deviations  can  be  ob¬ 
tained  as  a  function  of  azimuth  and  distance.  The  variation 
of  p  and  azimuth  deviations  shown  by  this  study  suggest  that 
the  station  deviations  depend  mainly  on  azimuth*  However,  it 
is  possible  that  they  also  depend  on  distance  to  some  extent. 
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D.  An  Analysis  of  a  Technique  for  the  Generation  of 

High  Resolution  Wavenumber  Spectra. 

The  technique  of  computing  high-resolution  wavenumber 
spectra  follows  the  method  used  by  Haney,  (1967),  The  basic 
method  consists  of  designing  a  multichannel  filter  whose  de¬ 
sired  output  is  an  impulse  at  time  t  =  t  and  spatial  lag  x  = 
V  where  tQ  is  usually  zero  and  xq  is  the  spatial  position  of 
one  of  the  seismometers.  The  method  consists  of  designing  a 
least-mean-square  error  "whitening"  filter  to  "whiten"  the 
data  in  f-k  space.  The  inverse  of  the  f-k  response  of  the 
filter  which  given  a  "white"  or  constant  f-k  spectrum  when 
applied  to  the  space-time  series  data  should  give  a  good  es¬ 
timate  of  the  true  f-k  spectrum  of  the  original  data.  This 
method  could  also  be  considered  as  equivalent  to  extending 
the  spatial  correlation  functions  (Burg,  1967). 

The  multichannel  whitening  filter  is  sensitive  to 
gain  inequalities.  One  method  of  compensation  is  to  norma¬ 
lize  the  spectral  matrix  in  a  manner  prescribed  by  Haney, 

(1967) . 

We  further  note  that  the  expression  derived  in  the 

i 

original  report  for  the  high-resolution  wavenumber  spectrum 
for  two  plane  waves  input  is  quite  similar  to  that  for  the 
spectrum  of  a  single  plane  wave,  with  the  exception  that 
P(f,k)  has  cross  terms,  and  N  is  replaced  by  R(  u )  in  the 
expansion  (c  +  N)  .  Some  of  the  results  for  a  single 
plane  wave  and  two  plane  waves  input  are  described  in  the 
next  section. 
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Figure  18.  The  WMO  array  response  to  an  infinite-velocity 
plane  wave. 


Figure  19  is  the  ordinary  f-k  spectrum  of  the  12  km/sec 
impulse  at  a  frequency  of  .3125  cps.  Figure  20  is  the  high- 
resolution  wavenumber  spectrum  (designed  on  the  center  seismo 
meter)  for  the  same  event  at  the  same  frequency.  Comparing 
Figure  20  with  Figure  19  we  see  that  the  high-resolution  wave- 
number  spectrum  at  .325  cps  indeed  shows  an  improvement  in 
resolution  over  the  ordinary  spectrum  (Figure  19).  Figure  20 
shows  a  signal  arriving  from  the  south  with  a  velocity  of 
slightly  greater  than  10  km/sec,  whereas  Figure  19  (the  ordi¬ 
nary  f— !..  spectrum)  shows  a  signal  coming  from  the  general 
direction  of  the  south  but  with  much  uncertainty  as  to  the 

exact  velocity  and  azimuth. 

Figure  21  is  the  high-resolution  wavenumber  spectrum 

based  on  the  output  of  a  WMO  sensor  at  a  frequency  of  .3125  cps. 

Comparing  Figure  21  with  Figure  20,  we  see  that  the  resolution 

of  the  high-resolution  f-k  spectrum  (Figure  21)  is  approximately 

equal  to  the  resolution  of  the  high-resolution  f-k  spectrum 

(Figure  20)  obtained  using  the  output  of  another  seismometer; 

however,  the  spectrum  in  Figure  21  is  "stretched"  along  the 

line  perpendicular  to  the  direction  of  the  line  connecting  the 

two  recording  positions.  Thus,  the  net  result  of  designing 

a  high-resolution  f-k  spectrum  on  a  seismometer  located  at 

x  =  x  is  to  introduce  a  cos(2tt  k»x  )  term,  which  causes  a 
~o  ~  . 

stretching  of  the  spectrum  along  a  direction  perpendicular  to 

to  the  vector  x  =  x  .  Figure  22  is  the  high-resolution  f-k 
~  ~o 

spectrum  averaged  over  all  of  the  seismometers  in  the  array, 
for  the  12  km/sec  event,  at  a  frequency  of  .3125  cps.  Com¬ 
paring  Figure  22  with  Figures  20  and  21,  we  see  that  in  the 
high-resolution  spectrum  which  has  been  averaged  over  all  of 
the  seismometers,  the  non-physically  meaningful  stretching  of 
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Figure  19.  The  ordinary  f-k  spectrum  of  a  theoretical 
12.0  km/sec  plane  wave  arriving  at  WMO  from 
the  south;  the  frequency  is  0.3125  cps. 
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Figure  21.  The  high-resolution  f-k  spectrum  for  the  12.0 
km/sec  plane  waves  at  frequency  at  0.3125  cps. 
The  reference  seismometer  is  located  at 
(.94,-. 17) ;  S/N  =  .03. 
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The  high-resolution  f-k  spectrum  for  the  12 
km/sec  plane  wave  at  .3125  cps  averaged  ove 
all  of  the  seismometers;  S/N  =  .03. 


the  spectrum  due  to  designing  the  high-resolution  spectrum 
on  an  off-center  seismometer  has  been  cancelled  out.  The 
spectrum  of  Figure  21  is  not  quite  as  peaked  as  that  of 
Figure  20,  although  the  difference  is  insignificant.  The 
spectrum  of  Figure  20  is  also  tilted  to  the  left  due  to  the 
asymmetry  of  the  WMO  array,  whereas  the  spectrum  of  Figure 
21  seems  to  be  slightly  "stretched"  in  a  northwardly  direction. 

Figure  23  is  the  ordinary  f-k  spectrum  of  LONGSHOT 
recorded  at  LASA  Fl  at  a  frequency  of  0.625  cps.  Figure  24  is 
the  "averaged"  high-resolution  f-k  spectra  at  the  same  fre¬ 
quency.  The  signal-to-noise  ratio  (c)  was  2.0.  At  0.C25  cps, 
both  the  averaged  high-resolution  (Figure  24)  and  the  ordinary 
spectrum  (Figure  23)  indicate  the  event  is  coming  from  the 
correct  direction,  but  the  velocity  is  somewhat  smaller  than 
expected  (11.6  km/sec,  vs.  14.1  km/sec,  from  the  J-B  traveltime 
tables) . 

Figure  25  is  a  beam-steered  (time-shift  and  sum)  power 
density  spectrum  of  LONGSHOT.  Nine  seismometers  in-line  in  the 
approximate  direction  of  the  event  were  time-shifted  and  sum¬ 
med.  The  power  spectrum  of  thu  sum  was  then  computed  and  dis¬ 
played  in  Figure  25.  This  demonstrates  that  there  is  enough 
power  at  2.5  cps  to  give  an  accurate  f-k  spectrum  at  that  fre¬ 
quency. 

Figure  26  is  the  array  response  of  thirteen  of  the 
seismometers  of  LASA  XE3.  Figure  27  is  the  ordinary  f-k 
spectrum  of  a  plane  wave  of  velocity  20  km/sec  and  amplitude 
1.0  plus  a  plane  wave  of  velocity  30  km/sec  and  amplitude 
0.2,  computed  at  0.625  cps.  Figure  28  is  the  averaged  high- 
resolution  f-k  spectrum  (S/N  =  2.0)  of  the  same  event  calculated 
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Figure  23.  Ordinary  f-k  spectrum  of  LONG  SHOT  at  LASA  sub¬ 
array  FI  at  .625  cps. 
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Figure  24.  High-resolution  f-k  spectrum  of  LONG  SHOT  at 
LASA,  FI,  averaged  over  all  st  ? 1  mometers  in 
the  array  at  .625  cps.  C=  2.i 
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Figure  26.  Array  response  of  the  partial  XE3  array  to  an 
infinite  velocity  plane  wave. 
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Figure  27.  Ordinary  f-k  spectrum  at  .625  cps  of  two  plane 
waves,  one  of  amplitude  1  at  20  km/sec,  North; 
and  one  of  amplitude  .2  at  30  km/sec,  South. 
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Figure  28.  High-resolution  wave number- spectrum  of  the  two 
plane  waves  at  .625  cps  averaged  over  all  of 
the  sensors  in  the  array  S/N  =  2.0. 


at  the  same  frequency.  Figure  29  is  the  high-resolution  f-k 
spectrum  (S/N  =  0.03)  of  the  same  event,  calculated  at  the 
same  frequency.  At  .625  cps  the  30  km/sec  event  would  not 
be  detected  by  either  the  ordinary  f-k  spectrum  (Figure  27) 
or  by  the  average  (S/N  =  2.0)  high-resolution  f-k  spectrum 
(Figure  28).  The  extreme  elongation  of  the  average  (S/N  = 
.03)  high-resolution  f-k  spectrum  (Figure  29)  might  make  one 
suspicious  about  the  presence  of  a  second  event,  but  the  con' 
elusion  would  depend  on  the  subjective  discrimination  of  the 
analyst. 

We  conclude  from  this  investigation  that  although 
the  high-resolution  frequency-wavenumber  spectrum  gives 
much  better  resolution  for  single  events,  the  technique  does 
not  impro /e  significantly  the  capability  of  a  seismic  array 
to  detect  multiple  time-overlapping  events  from  different 
azimuths.  A  possible  method  of  detecting  a  small  event  in 
the  presence  of  a  large  event  is  to  design  an  approximate 
pre-whitening  filter  in  the  time  domain  to  approximately 
"whiten"  the  f-k  spectrum  of  the  data,  convolve  the  filter 
with  the  time  series  data,  and  then  proceed  with  either  or¬ 
dinary  or  high-resolution  f-k  analysis.  The  fact  that  the 
pre-whitening  filter  need  be  only  an  approximate  whitening 
filter  facilitates  its  design. 
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Figure  29.  High-resolution  wavenumber  spectrum  of  the  two 
plane  waves  at  .625  cps  averaged  over  all  of 
the  sensors  in  the  array.  S/N  =  .03. 


Haney,  W.P.,  1967,  Research  on  high-resolution  frequency- 
wavenumber  spectra,  Special  Scientific  Report  No.  2: 

Dallas,  Texas  Instruments,  Inc.,  Contract  AF33  (557) -16678. 

McCowan,  D.W.,  1966,  Finite  Fourier  transform  theory  and 
its  application  to  the  computation  of  convolutions,  cor¬ 
relations,  and  spectra,  Seismic  Data  Laboratory  Report  No. 

168  Revised:  Alexandria,  Virginia,  Earth  Sciences,  A  Teledyne 
Company. 

McCowan,  D.W.  and  Lintz,  P.R.,  1968,  High-resolution  fre¬ 
quency-wavenumber  spectra,  Seismic  Data  Laboratory  Report 
No.  206,  Alexandria,  Virginia,  Earth  Sciences,  A  Teledyne 
Company. 

E.  Vertical  Array  Frequency-Wavenumber  (v-f-k)  Spectra 

of  Synthesized  Rayleigh  Wave  and  P-Wave  Noise. 

In  recent  reports,  Sax  (1967)  has  discussed  v-f-k 
(vertical  array  frequency-wave  number)  spectra  of  the  noise 
fields  preceding,  underlying  and  following  high  frequency 
teleseisniic  P  waves.  In  some  cases,  difficulties  in  relating 
v-f-k  spectra  of  seismic  data  to  wavemodes  were  encountered. 

Model  studies  "\re  therefore  undertaken  in  this  study  to  inves¬ 
tigate  v-f-k  p  it ter ns  of  synthesized  Rayleigh  wave  and  body 
wave  noise.  The  purpose  is  to  identify  v-f-k  patterns  of  syn¬ 
thesized  Rayleigh  wave  and  body  wave  noise.  The  purpose  is  to 
identify  v-f-k  patterns  associated  with  random  Rayleigh  waves 
or  body  wave  components . 

McCowan  (1968)  described  the  method  for  computing 
v-f-k  spectra.  Harkrider's  (1964)  computer  program  was  used 
to  compute  the  Rayleigh  wave  modes  described  in  the  following 
results. 

Figure  34  shows  the  geological  structure  used  to 
compute  the  Rayleigh  wave  modes.  The  layering  is  character¬ 
ized  grossly  by  a  halfspace  velocity  of  6  km/sec  and  inter- 
bedded  5.0  km/sec  granite  between  1.5  and  2.0  km.  Only  the 
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trapped  modes  produced  by  the  interbedded  low  velocity  zone 
are  considered.  Shear  velocities  were  estimated  by  assuming 
a  constant  Poisson's  ratio  of  .27. 

Amplitude-depth  curves  of  the  vertical  component  of 
displacement  are  shown  in  Figure  30  for  the  fundamental  mode; 
Figure  31,  first  higher  mode;  Figure  32,  second  higher  mode; 
and  Figure  33  the  third  higher  mode.  The  low  frequency  limit 
of  each  propagation  mode  is  determined  by  half-space  radiation 
of  modes  with  phase  velocities  greater  than  the  shear  velocity 
of  a  granite  half-space. 

The  noise  which  was  subjected  to  v-f-k  analysis  was 
synthesized  by  feeding  random  numbers  into  a  set  c  v  constant- 
Q  filters  of  .25  cps  effective  bandwidth.  In  order  to  cover 
a  broad  band  of  frequencies,  the  center  frequency  of  each 
filter  was  spaced  .25  cps  apart  starting  at  .25  cps  and 
ending  at  3.0  cps.  An  independent  realization  of  random 
numbers  was  fed  into  each  filter.  The  resultant  broad  band 
noise  was  the  sum  of  the  random  uncorrelated  outputs  of  each 
filter.  The  noise  at  the  surface  of  the  earth  was  normalized 
to  be  white,  i.e.,  with  a  constant  power  spectrum.  The  simu¬ 
lated  noise  at  depth,  z ,  was  multiplied  by  a  coefficient 

corresponding  to  the  n'th  propagation  mode  H  (z.,f  )  at  each 

n  l  m 

^re<?uency*  The  band— limited  wave  components  were  summed 

over  all  the  center  frequencies,  f  ,  to  obtain  at  each  level, 
zi<  a  broad  band  multichannel  realization  of  noise  in  the  n’th 
propagation  mode.  The  noise  is  the  convolution  of  a  white 
random  source  N‘(t,f  )  with  the  response  of  a  narrow  band 

F(t,f^)  which  peaks  at  f^„  These  narrowband  realizations 
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AMPLITUDE  (NORMALIZED) 


AMPLITUDE  (NORMALIZED) 


Figure  32.  Frequency- Amplitude-Depth  Relations 
For  Higher  Rayleigh  Modes  at  Station 
APOK  Second  Higher  Mode 


AMPLITUDE  (NORMALIZED) 


Figure  33.  Frequency-Amplitude- Depth  Relations 
For  Higher  Rayleigh  Modes  at  Station 
APOK  Third  Higher  Mode 
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Figure  34.  Velocity  and  Density  Sections  at  Site 
APOK  (Apache,  Oklahoma) 


m 


are  summed  by  weighting  with  H  (z.,f  ),  which  for  Rayleigh 

n  1  m 

waves  is  previously  cited.  Similarly,  body  wave  noise  was 

synthesized  by  using  the  halfspace  eigenfunction  cos 

(2it  ,  where  a  is  the  halfspace  velocity  of  P-waves. 

Figures  35-38  show  the  v-f-k  spectra,  for  k  >  o,  for 
the  four  Rayleigh  wave  modes  depicted  in  Figures  30  to  33. 

These  are  shown  as  if  measured  by  a  typical  six-element 
vertical  array  approximately  3  km  long,  with  each  channel 
filtered  by  an  operator  which  prewhitens  the  noise  at  the 
surface.  The  synthetic  data  were  generated  by  means  of  a 
random  number  generating  algorithm  and  narrow  band  filters 
as  previously  described.,  The  synthetic  data,  if  compared  to 
observed  data,  are  subject  to  the  same  array  response,  thereby 
eliminating  some  of  the  problems  of  array  response  implicit  in 
the  identification  of  waves.  This  follows  as  a  consequence  of 
using  the  same  array  geometry,  sample  size  and  rate,  smoothing, 
etc.,  for  both  sets  of  data. 

The  amplitude  spectra  of  the  higher  modes  (Figures 
36-38)  are  dominated  by  peaks  at  non-zero  wave  numbers  run¬ 
ning  parallel  to  the  frequency  axis.  This  reflects  a  fixed 
number  of  zero  crossings  as  a  taper  which  splits  the  wave 
number  response  into  two  nearly  parallel  peaks  similar  to  the 
way  a  cosine  taper  splits  a  dc  time  response  into  two  discrete 
spectral  peaks  at  the  frequency  of  the  taper.  Compare  the 
first,  second,  and  third  higher  modes:  the  first  higher  mode  is 
distinguishable  on  the  basis  of  wave  number  splitting,  but  the 
second  and  third  are  indistinguishable  from  each  other  on  this 
basis,  having  very  similar  v-f-k  patterns.  After  comparing  the 
displacement  curves  (Figures  32  and  33)  for  the  second  and 
third  higher  mode,  we  note  that  any  apparent  difference  in  wave 
number  response  is  below  the  resolution  of  the  v-f-k  processor. 
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VFK  SPECTRUM 


Figure  35.  Fundamental  Rayleigh  Mode  -  Vertical 
Component  -  Station  APOK 
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Figure  37. 


2nd  Higher  Rayleigh  Mode  - 
Component  -  Station  APOK 


Vertical 
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Figure  39  shows  vertical  P  wave  noise.  The  noise  is 
generated  by  using  the  halfspace  vertical  displacement,  cos 

2  it  f  x/a,  where  a  is  an  average  half  space  compressional 
velocity.  Figure  40  is  the  v-f-k  spectrum  corresponding  to  sn 
up-  and  down-tra veiling  impulse  with  the  same  average  half¬ 
space  velocity  used  for  Figure  39.  In  both  cases,  the  shaded 

3  db  peaks  indicate  an  apparent  vertical  phase  velocity  (c  = 
f/k  =  a)  corresponding  to  the  average  half-space  velocity,  a. 
The  array  responses  differ  considerably  between  the  P  impulse 
model  of  Figure  40  and  the  long  random  P  wave  function  of 
Figure  39.  The  time  series  resulting  in  Figure  40  would  be 
immediately  obvious  as  a  P  impulse.  The  random  time  series 
resulting  in  Figure  39  looks  almost  totally  incoherent  between 
channels.  It  cannot  be  recognized  as  P  noise  unless  we  refer 
to  Figure  39,  in  which  case  it  is  obvious. 

Comparing  the  v-f-k  analysis  of  higher  modes  on 
Figures  37  and  38  with  the  P  noise  on  Figure  39,  we  see  that 
P  noise  above  1.8  cps  is  barely  distinguishable  from  the  second 
and  third  higher  mode  Rayleigh  waves.  It  is  reasonable  to 
anticipate  some  difficulty  in  designing  filters  to  separate 
a  narrow  band-limited  P  signal  from  noise  containing  higher¬ 
mode  Rayleigh  waves.  This  result  is  significant,  as  it  has 
been  amply  demonstrated  that  higher  modes  and  possibly  leaky 
modes  and  P  waves  are  significant  components  of  the  ambient 
noise,  especially  above  .5  cps. 

The  results  of  the  study  indicate  that  theoretical 
v-f-k  analysis  of  wave  modes  should  be  extended  to  include 
deep  crustal  trapped  modes  which  are  almost  certainly  impor¬ 
tant  components  of  ambient  seismic  noise.  Further,  high- 
resolution  v-f-k  analysis,  McCowan  and  Lintz  (1968),  plus 
the  aforementioned  extended  array  concept  could  improve  the 
array  response  due  to  sampling  limitations  imposed  by  the 
array  geometry. 
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f  (cps) 


0.15  0.3  0.45  0.6  0.75 

k  (c/km) 

VFK  SPECTRUM 


40.  Teleseismic  P-wave  Signal  (white) - 
Vertical  Component  -  Station  APOK 


Table  IX  summarizes  part  of  the  seismic  data  compiled  from 
these  events . 


Other  projects  completed  included  photo  albums  for 
two  events  and  travel-time  computations  for  NTS  and  other 
sites.  In  addition  complexities  were  computed  for  a  number 
of  special  events. 


Name  Event  Date  Medium 
SCROLL  23  April  1968  Vitric  Tuff 


No. of  stations  recording 
Magnitude  SP  Signals  LP  Signals 


4.45  +  0.44 


BOXCAR  26  April  1968  Rhyolite  6.42  +  0.45  23  22 


TABLE  IX.  Seismic  Data  from  Two  Nuclear  Explosions 


III.  SUPPORT  AND  SERVICE  TASKS 

A.  VELA-Uniform  Data  Services 

As  part  of  the  contract  work  statement,  the  SDL  pro¬ 
vided  one  or  more  of  the  following  support  and  service  functions 
for  VSC  and  other  VELA  participants: 

-  copies  of  16  and  35  mm  film. 

-  playouts  of  earthquakes  and  special  events. 

-  copies  of  composite  analog  tapes. 

-  composite  analog  tapes  of  special  events. 

-  use  of  1604  computer  for  checking  out  new  programs  or 
running  production  programs. 

-  copies  of  digital  programs. 

-  digitized  data  in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604. 

-  running  SDL  production  programs,  such  as  power  spectral 
density  and  array  processing  on  specified  data. 
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-  digital  x-y  plots  of  power  spectra  or  digitized  data. 

-  signal  reproduction  booklets. 

-  space  for  visiting  scientists  utilizing  SDL  facilities 
to  study  data  and  exchange  information  with  SDL 
personnel. 

During  this  report  period,  56  such  projects  were  com¬ 
pleted  and  the  9  organizations  receiving  these  services  are 
listed  in  Appendix  A. 

B.  Data  Library 

The  data  library  contains  approximately  8130  digitized 
seismograms,  237  digital  computer  programs  and  298  composite 
analog  magnetic  tapes,  all  available  for  use  by  the  VELA-Uniform 
program. 

The  following  additions  were  made  during  this  period: 

1.  Digital  Seismograms  -  730  including 

-  data  from  two  explosions 

-  one  earthquake  recorded  at  various  stations 

2.  LASA  Data  -  12  tapes 

-  there  are  a  total  of  1408  digital  tapes  in  the  library 
including  1004  field  tapes.  There  is  also  a  master 
calibration  tape  which  contains  the  magnification 
(digital  counts  per  millimicron)  of  each  sensor  for 
every  subarray.  These  magnifications  have  been  com¬ 
puted  for  all  calibration  tapes  currently  in  house. 

As  each  new  calibration  is  received,  it  is  routinely 
run  through  the  new  program  CALIBR  and  added  to  the 
master  tape. 

3.  Digital  Programs  -  12  including 

5CRTAPE  -  To  make  and  update  BCD  tapes  containing 

standard  production  programs  and  subroutines. 


* 


-30- 


COOLCEP  -  This  program  computes  and  displays  power 

spectra  and  cepstra.  A  plot  tape  is  generated 
with  explanatory  labels.  Provision  has  been 
made  for  tapering  and  high-pass  filtering 
(littering)  the  power  spectrum  prior  to  com¬ 
puting  the  cepstrum.  The  Cooley-Tukey  method 
is  used  to  obtain  high  speed.  The  program  is 
designed  for  multiple  runs,  reading  data  cards 
until  and  end-of-file  is  reached  on  the  card 
reader. 

BMSTR  -  Beamforms  the  data  channels  of  an  input  seismogram 
and  computes  the  power  spectra  from  selected  parts 
of  the  resulting  trace. 

BMSTEER  -  Generates  a  set  of  MCF  beam-steer  deltas  and 

stores  them  on  magnetic  tape  in  a  format  suit¬ 
able  for  processing  by  program  M2 20  MCFONPT 
(revised) .  The  time  shifts  are  calculated 
geometrically  from  the  array  coordinates  which 
are  stored  on  tape  and  a  velocity  and  back 
azimuth  read  in  from  cards. 

MXTFOCAL  -  Computes  magnification  levels  in  counts/microns 
from  the  TFO  calibration  data  on  the  original 
TFO  acquisition  system  tape. 

CALSETUP  -  This  subroutine  removes  the  mean  and/or  linear 
trend,  finds  the  calibration  and  computes  mag¬ 
nification  level  in  counts/microns. 

RESTRA  -  Determine  the  coefficients  of  the  linear  equation 
y  =  a  +  Bx  or  alternately  Log  (Z)  =  c  +  d  Log  (x) . 
The  user  has  the  option  to  choose  either  equation. 
Moreover,  the  user  can  decide  if  correction 
factors  are  to  be  applied  to  the  y's  or  the  Log 
(Z)'s. 

REGRESS  -  This  program  generates  station  correction  factors 
which  are  used  as  input  to  program  RESTRA. 
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PREFILT  -  This  program  writes  a  save  tape  of  a  MCF  pre¬ 
filter  suitable  for  processing  by  program  M220 
MCFONPT  (revised) .  The  prefilter  coefficients 
are  read  from  cards, 

MULMATCH  -  Performs  multichannel  matched  filtering  on 

array  data.  The  program  produces  printed  and 
plotted  outputs  displaying  the  filter,  the 
individual  filtered  traces,  and  the  sum. 

NORSARLP  -  Processes  a  Norsar  tape  (merged  by  Lincoln 
Laboratories)  containing  multiplexed  long- 
period  seismic  data,  and  forms  an  unpacked 
SDL  library  tape. 

GIVHOUSE  -  This  Fortran-63  subroutine  calculates  the  eigen¬ 
values  and  eigenvectors  of  a  real,  symmetric 
matrix.  The  Givens-Householder  method  is  used 
to  reduce  the  matrix  to  a  tri-diagonal  form. 

The  method  of  bisection  is  then  used  to  obtain 
the  eigenvalues  and  inverse  iteration  to  obtain 
eigenvectors.  The  eigenvectors  are  normalized 
to  have  unit  magnitude.  An  attempt  is  made  to 
provide  linearly-independent,  but  not  necessarily 
orthonormal,  eigenvectors  for  degenerate  eigen¬ 
values.  Provision  also  exists  to  calculate  less 
than  the  normal  number  of  eigenvalues  or  eigen¬ 
vectors.  The  eigenvalues  are  returned  in  order 
of  decreasing  magnitude. 

Analog  Composite  Tapes  -  5  including 

a.  Made  by  SDL 

STUTZ 

REX 

RED  HOT 

b.  Made  by  Geo tech 

KNOX 

STINGER 
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C.  Data  Compression 

This  is  a  continuing  routine  operation,  and  production 
is  maintained  at  the  level  needed  to  meet  the  requirements  of 
the  field  operation  (LRSM  and  U.S.  Observatories)  and  the 
Seismic  Data  Laboratory.  For  this  period  2100  tapes  were 
compressed. 

D.  Automated  Bulletin  Process 

April,  May  and  June  1968  LRSM  and  Observatory  Bulle¬ 
tins  were  processed  during  this  report  period  and  forwarded  to 
Geotech,  a  Teledyne  Company,  for  checking  and  publication. 
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April  -  June  1968 


nnnn^aHor.g  Receiving  SDL  Data  Services 


Vitro 

General  Atronics  Corporation 
Geotech,  A  Teledyne  Company 
Oregon  State  University 
Texas  A&M  University 
S.W.  Center  for  Advanced  Studies 
U.S.  Coast  St  Geodetic  Survey 
University  of  North  Carolina 
St.  Louis  University 
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